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'  Hydrostatic  pressure,  compression  under  different  confining  pressures, 
proportional  loading  and  uni  axial -strain  tests  were  conducted  on  a  dry  tuff 
from  the  Nevada  Test  Site  and  dry  Kayenta  sandstone  from  the  Nixed  Company 
Site  near  Grand  Junction,  Colorado.  Both  rocks  exhibit  behavior  which  Is 
characteristic  of  very  porous  geological  materials;  namely.  Increasing  shear 
strength  as  a  function  of  confining  pressure,  large  volume  compaction 
followed  by  volume  dilatancy  under  selected  loadings  with  enhanced  compact  on, 
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over  some  pressure  ranges,  in  the  presence  of  shear  deformation.  Both 
materials  also  exhibit  strain  work-hardening.  A  “cap  model"  is  applied 
to  represent  the  mechanical  behavior  of  these  two  materials.  This  model 
accounts  for  all  of  the  aforementioned  characteristics  of  both  the  tuff 
and  the  sandstone,  and  provides  a  means  to  represent  the  path  dependence 
of  their  mechanical  properties. 
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I.  INTRODUCTION 


A  shallow-buried  structure,  such  as  an  ABM  silo,  may  be  made 
relatively  safe  from  nuclear  and  thermal  radiation  and  from 

the  direct  effects  of  air-blast.  The  primary  vulnerability,  therefore, 
is  associated  with  the  motion  of  the  surrounding  earth.  In 
addition  to  the  intense  direct  shock  in  the  ground  that  is  responsible 
for  the  crater  formation  from  a  surface  burst,  ground  motions  are 
induced  by  the  passage  of  an  air  blast  wave  ’long  the  surface.  For 
most  surface  or  shallow  buried  structures,  this  air-i'duced  ground 
shock  is  of  great  significance  since  it  is  extended  to  large  distances, 
while  the  direct  ground  shock  is  more  rapidly  attenuated  in  parsing 
through  the  intervening  earth  mass. 

In  the  past,  direct  ground  shock  was  calculated  using  hydro- 
dynamic  models  while  air-induced  ground  shock  effects  were  approximated 
by  means  of  linear  elastic  models.  Experimental  data  has  shown  that 
these  models  do  not  represent  geologic  materials  well,  and  hence 
more  realistic  models  are  required.  To  calculate  both  the  intense 
direct  ground  shock  and  the  air-induced  ground  shock  requires  a 
description  of  the  material  behavior  in  the  stress  range  extending 
from  Megabar  stresses,  which  occur  near  ground  zero,  down  to  very  low 
stresses,  which  occur  at  large  distances  from  ground  zero.  It  is 
probably  the  stress  region  between  the  very  high  stresses,  which  may 
be  reasonably  well  represented  by  a  hydrodynamic  model,  and  the  extremely 
low  stresses  which  may  reasonably  well  be  represented  by  linear  elasticity, 
that  require  the  most  complicated  constitutive  model.  In  particular,  the 
stress  region  from  the  linear  elastic  region  to  ten's  of  kilobars  requires 
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the  most  attention.  This  is  the  stress  region  covered  by  the  work 
here. 

Stress-strain  data  for  two  geological  materials  have  been 
obtained  from  a  variety  of  load-path  tests.  The  two  geologic  mater¬ 
ials  were  selected  to  represent  two  extremes  in  very  porous  geologic 
materials;  one  with  relatively  "weak  cementing  of  grains"  and  the 
other  with  relatively  "strong  cementing  of  grains."  The  tests  serve 
two  major  purposes:  (i )  to  give  the  general  characteristic  stress- 
strain  behavior  in  order  to  build  mathematical  models  that  qualita¬ 
tively  fit  the  behavior  at  the  appropriate  stress  levels,  and  (ii) 
to  provide  quantitative  data  from  which  the  required  material 
constants  in  the  mathematical  models  can  be  determined. 

The  observed  data  are  fit  with  a  "cap  model."  This  elastic- 
plastic  model  accounts  for  the  general  features  and  fits  the  hydro¬ 
static  pressure  and  compression  response  under  different  confining 
pressures  very  well.  Some  shortcomings  of  the  "cap  model"  as 
applied  here  are  discussed  and  some  in-situ,  real  rock  mass  effects 
are  mentioned. 
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II.  EXPERIMENTAL  RESULTS 


Experimental  Techniques 
Loading  Presses 

Laboratory  tests  were  conducted  on  smai?  >olid  cylindrical 
samples  using  servocontrolled  testing  machines.*  The  samples  were 

f 

subjected  to  fluid  confining  pressure  and  to  axial  compression 

simultaneously.  By  programing  the  application  of  the  load  and  the 

application  of  the  fluid  pressure  independently,  a  variety  of  load- 
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path  tests  were  conducted.  Figure  la  shows  the  servocontrolled 
loading  presses  used  for  tests  here. 

Sample  Size 

Test  specimens  were  made  by  diamond  coring  and  finish  grinding 

samples  to  final  dimensions  of  1.9  cm  (.75  inch)  diameter  by  3.8  cm 

(1.50  inch)  long.  The  specimens  were  covered  with  .25  mm  thick 

urethane  sheet  and  metal  end  caps  were  epoxied  to  the  specimen  ends. 

The  urethane  sheet  prevented  fluid  intrusion,  but  carried  an  insigni- 
2 

ficant  load.  More  details  of  test  specimen  preparation  are  contained 
in  Appendix  B. 

Measurement  of  Stresses  and  Strains 

All  stresses  and  strains  were  measured  inside  the  pressure 
1  2 

vessel,  *  and  therefore  corrections  for  seal -friction,  alignment, 
pressure  and  temperature  effects  were  either  not  needed  or  were 
measured  directly.  Pressure  was  measured  with  a  manganin  pressure 
coil;  load  was  measured  using  a  load-cell  adjacent  (axially  in-line 
with)  to  the  test  specimen.  Axial  and  transverse  strains  were 
measured  using  cantilevered  displacement  measuring  devices,  precalibrated 
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Figure  lb.  Test  specimen  with  cantilever 
Figure  1.  Test  Apparatus. 


inside  the  pressure  vessel.  The  axial -cantilever  displacement  sensor 
provided  a  measure  of  average  specimen  shortening,  while  the  trans¬ 
verse-cantilever  sensor  provided  measurement  of  average  diameter 
change  at  two  diameters.  Figure  lb  shows  the  test  specimen  with  the 
axial  and  transverse  strain  sensors  attached. 

Aoouracu  of  Measurements 

The  pressure  measurement  is  accurate  to  ±10  bars  and  load  to  ±10 
bars.  The  cantilever  displacement  sensors  give  cantilever  displacements 
accurate  to  ±5  x  10  4  mm.  Average  specimen  strains  are  measured  accur¬ 
ate  to  better  than  500  microstrains  (5  x  10’4  strain).  Uncertainty 
in  measuring  strains  is  due  to  the  specimen  jacket  materials  and  the 
specimen -jacket  interface.  Repeatability  of  several  tests  is  shown 
by  scatter  bars  on  the  figures. 

All  tests  were  conducted  at  room  temperature;  the  loading  was 
such  as  to  give  strain  rates  of  about  10"4 /second.  Details  of  the 
techniques  are  contained  in  Appendix  A.  Throughout  this  work  the 
volume  strain  is  calculated  according  to  the  small  strain  approximation, 
i.e.,  A V/VQ  =  ej  +  2e3  ,  where  ei  and  £3  are  the  axial  and  transverse 
strains,  respectively.  The  error  involved  in  this  approximation  is 
shown  in  Appendix  A  to  be  less  than  the  scatter  of  the  individual 
test  data.  This  approximation  is  also  commensurate  with  mathematical 
models  for  porous  rock  behavior,  such  as  the  one  given  in  Section  IV. 

Material  Description 

The  two  materials  tested  were  an  ash-fall  tuff  from  the  Nevada 
3 

Test  Site  and  a  competent  sandstone  from  the  Mixer*  Company  high- 
energy  explosive  site  near  Grand  Junction,  Colorado.4  Physical 
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properties  are  given  below;  more  detailed  descriptions  of  the  rocks 
are  included  in  Appendix  B.  All  tests  were  conducted  on  "dry"  samples 
The  effect  of  water  content  or  the  mechanical  properties  of  tuff  and 
sandstone  is  beyond  the  scope  of  this  report  and  remains  to  be  fully 
investigated. 


TABLE  1.  PHYSICAL  PROPERTIES  OF  ROCKS  TESTED 


Dry  Density 
ga/cm3 

Grain  Density 
gm/an3 

Total  Porosity 
t  of  Volume 

Grain  Size 

Pore  Size 

Tuff 

1.7 

2.46 

31 

10"7mm 

(fragments  to 

to  >  1  Ml 

5mm) 

Sandstone 

2.0 

2.66 

25 

0.2  -  0.3  am 

0.1  mi 

Experimental  Results  on  Tuff 

Loading  Paths 

A  number  of  different  loading  paths  were  selected,  including: 
(1)  hydrostatic  pressure,  (2)  compression  under  different  constant 
confining  pressures  from  zero  to  8  kbars  ( tri axial* stress ) ,  (3) 
uniaxial-strain  and  (4)  proportional -stress.  During  these  tests  the 
stress  state  Is  given  below  in  Table  2,  where  the  1  direction  refers 
to  the  axis  of  the  cylindrical  specimen  and  2  and  3  are  transverse 
orthogonal  directions.  Compressive  loading  and  shortening  strains 
are  taken  as  positive.  All  stress-strain  values  are  engineering 
stress-strain,  without  correcting  for  geometry  changes. 


TABLE  2.  STRESS  STATES  FOR  VARIOUS  LOADING  PATHS 


Stress -State 

Hydrostatic 

0|  ■  02  ■  OJ 

Trl axial  Stress 

0|  >  O;  •  0} 

Uniaxial  Strain 

lax  >  o2  ■  o3\ 

\e2  *  *3*0/ 

Proportional  Stress 

/ax  >  o2  -  o3  \ 

\oi/o3  ■  const./ 

The  loading  paths  are  shown  In  Figure  2,  where  the  load,  i.e., 
stress  difference  (01  -  o3),  has  been  plotted  versus  the  mean  normal 
stress  for  each  of  the  different  kinds  of  tests.  Each  line  repre¬ 
sents  the  average  of  three  to  five  tests,  with  the  loading  preceding 


Figure  2.  Loading  paths  for  dry  NTS  ash-fall  tuff. 
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to  the  end  of  the  line  as  shown.  Unloading  follows  identically  the 

loading  paths  except  for  the  uniaxial- strain  tests.  The  hydrostatic- 

pressure  tests  lie  along  the  mean  normal  stress  axis,  while  triaxial- 

stress,  uniaxial -strain  and  proportional -stress  tests  are  shown  with 

different  symbols.  The  "shear-stress  limit"  curve  is  an  estimate  of 

what  is  conventionally  called  the  "failure  envelope."  This  limit- 

curve  constitutes  a  maximum  shear-stress  (load  here)  that  can  be 

5 

reached  for  any  path  of  loading;  i.e.,  stress  states  beyond  this  curve 
are  not  possible.  The  locus  of  this  curve  is  determined  by  the 
asymptotic  limit  to  the  stress  difference  oj  -  o3  as  a  function  of  the 
axial  and  transverse  strains  in  a  tri axial  compression  test. 

Hyr?o8tatic  Teste 

The  hydrostat  for  this  tuff  is  snown  in  Figure  3,  where  the 
volumestrain,  aV/V0,  is  obtained  as  the  measured  axial  strain  plus 
the  two  transverse  strains.  Axial  and  transverse  strains  were  measured 
and  the  results  suggest  isotropic  response.  Even  though  large  volume 
strains  occur,  approximately  ten  percent  porosity  still  exists  at  8  kbars 
pressure.  Apparently  much  higher  pressure  is  required  to  completely  compact 
all  porosity.  Several  unloadings  and  reloadings  are  shown.  Considerable 
recovery  occurs  as  is  noted  by  the  bend-back  of  the  hydrostat  when  un¬ 
loading  from  8  kbars,  for  example. 

Triaxial-Stresa  Teste 

The  stress-strain  curves  for  compression  at  different  confining 
pressures  are  shown  in  Figure  4.  The  values  of  the  confining 
pressure,  o3,  are  shown  in  the  figure.  The  axial  shortening  strain 
and  the  transverse  extension  strain  during  the  compressive  loading 


are  shown.  From  this  figure,  shear  strain  and  volume  strain  can  be 
obtained  (and  are  shown  later),  as  well  as  the  apparent  elastic 
constants,  Young's  modulus  and  Poisson  ratio.  A  large  increase  in 
Young's  modulus  is  shown  with  increasing  confining  pressure.  A 
transition  from  brittle  behavior,  for  the  unconfined  compression 
test,  to  ductile  behavior  for  all  the  confined  tests,  occurs  when 
brittle  refers  to  an  abrupt  load  drop  after  some  maximum  stress 
and  ductile  refers  to  the  ability  to  undergo  large  crushup  strains 
without  a  load  drop.  Work-hardening,  or  strain  strengthening,  occurs 
to  very  large  strains  for  the  confined  tests.  Unloading  occurs  along 
a  much  steeper  modulus  than  loading. 

Uniaxial-Strain  and  Proportional-Stre88  Teste 
The  stress-strain  curves  for  uniaxial -strain  and  two  different 
proportional  stress  paths*  oj/o 3  *  1.88  and  oi/o3  =  2.81,  are  shown 
in  Figure  5.  Here  the  axial  and  transverse  strains  are  the  "total" 
strains  starting  from  zero  stress,  and  from  these  the  shear  and 
volume  strains  can  again  be  obtained. 

The  oj/03  *  1.88  proportional -stress  and  the  uniaxial-strain 
tests  follow  relatively  similar  paths,  while  the  0l/o3  =  2.81  propor¬ 
tional-stress  tests  show  much  less  work  hardening  and  greater 
transverse  strain.  Considerable  recovery  of  the  axial  strain  occurs 
for  the  proportional -stress  load  (unload)  paths,  similar  to  the  hydro- 
stat.  Figure  3.  The  uniaxial -strain  path  shows  somewhat  less  recovery 
on  unloading  than  the  triaxicl  or  hydrostatic  tests. 

*  These  two  stress  paths  are  chosen  by  drawing  the  approximate  location 
of  the  desired  stress  path  on  graph  paper  and  then  accurately  measur¬ 
ing  the  slope  followed  after  the  test  is  completed.  There  is  no 
special  significance  to  the  particular  values  of  used  here. 
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Figure  5.  Uniaxial-strain  and  proportional -stress 
tests  on  dry  NTS  ash-fall  tuff. 


Experimental  Results  on  Sandstone 

Loadvig  Paths 

The  different  load  paths  selected  were  similar  to  those  for  the 
tuff.  Figure  6  shows  the  load  versus  mean  normal  stress  for  the 
different  paths;  each  path  represents  the  average  of  three  to  five  tests. 
A^ain  the  hydrostatic  pressure  loading  lies  along  the  mean  normal  stress 
axis,  while  other  paths  are  shown  with  different  symbols.  The  sandstone 
shows  a  shear-stress  limit  relatively  similar  to  the  tuff,  up  to  mean 
normal  stress  of  about  4  kbars.  For  higher  mean  normal  stresses  the 
shear  strength  appears  to  be  plateauing,  increasing  only  slightly  with 
tsean  normal  stress.  The  uniaxial -strain  path  shown  in  Figure  6  for  sand¬ 
stone  differs  significantly  from  the  nearly  linear  response  shown  in 
Figure  2  for  the  tuff. 
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Hydrostatic  Teste_ 

The  hydrostat  for  the  dry  sandstone  Is  shown  In  Figure  7.  As 
for  tne  tuff,  the  axial  and  transverse  strains  were  approximately 
equal,  suggesting  isotropic  response  The  sandstone,  unlike  the 
tuff,  shows  a  “foot."  or  elastic  compaction,  up  to  about  .3  to  .4 
kbars  pressure.  This  is  believed  due  to  closing  of  micro-cracks 
which  exist  in  the  sandstone.  Large  permanent  compaction  occurs 
for  pressures  above  about  4  kbars.  However ,  the  volume  compacts 
of  about  15*  at  8  kbars  pressure  Indicates  that  not  all  the 
Initial  25*  void  space  has  been  compacted,  even  at  this  high 
pressure.  The  competent  Intact  sandstone  specimen  becomes  "crushed" 
broken  for  pressures  above  about  4  kbars  and  appears  to  be  loose 


"sand"  after  unloading. 


Figure  7.  Hydrostat  for  dry  sandstone. 
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Triaxial-Stress  Tests 


The  stress-strain  curves  for  compression  at  different  confining 
pressures  are  shown  in  Figure  8,  where  axial  shortening  strain  and 
transverse  extension  strain  are  plotted  versus  the  compressive  load. 
The  confining  pressures  are  shown  on  each  curve  as  o3 .  Increase  in 
the  Young's  modulus  is  shown  with  increasing  confining  pressure  and 
a  transition  from  brittle  behavior  for  the  unconfined  test  to  ductile 
behavior  for  all  the  confined  tests  occurs.  Work  hardening,  or 
strain  strengthening,  is  shown  (for  the  confined  tests)  to  very 
large  shortening  strains.  Unloading  occurs  along  a  much  steeper 
modulus  than  loading. 


TRANSVERSE  AXIAL 

STRAIN, % 

Figure  8.  Tri axial  compression  tests  for  dry 
sandstone. 
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ihnj./iarLal-Stmin  and  Proportional-Stress  Tests_ 

The  stress-strain  curves  for  uniaxial -strain  and  proportional  - 
stress,  0l/o3  »  1.88  and  0l/o3  =  3.18,  loadings  are  shown  in 
Figure  9.  Axial  and  transverse  strains  are  "total"  strains  starting 
from  zero  stress.  From  these  curves,  shear-strain  and  volume- 
strain  can  be  determined  and  are  presented  later. 

Just  as  for  the  tuff,  the  a./o,  =  1.88  proportional -stress  test 
and  the  uni  axial -strain  tests  follow  relatively  similar  paths.  The 
higher  stress  ratio  proportional -stress  test  shows  less  (no)  work 
hardening  and  larger  transverse  strain. 


Figure  9. 


Uniaxial -strain  and  proportional 
stress  tests  on  dry  sandstone. 
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III.  GENERAL  FEATURES 


Hydrostatic  Compression 

Hydrostatic  compaction  at  low  pressure,  jp  to  hundreds  of  bars 
for  the  dry  tuff  and  about  4  kbars  for  the  dry  sandstone,  will  give 
elastic,  but  generally  non-linear,  response  with  no  or  little 
measurable  permanent  set.**  The  initial  bulk  modulus  in  this  approxi¬ 
mately  elastic  region  would  agree  with  that  calculated  from  ultrasonic 
velocity  measurements  of  shear  and  longitudinal  velocities, 
except  that  (i)  strain-rate  effects  are  neglected  and  (ii)  ultrasonic 
pulses  are  of  a  micro-strain  magnitude  whereas  even  for  pressures  of 
bars  the  strains  are  orders-of-magnitude  greater  and  hence,  may  result 
in  different  deformation  mechanisms.  The  calculated  elastic  constants 
from  measurement  of  shear  and  longitudinal  velocities7  (of  dry  sand¬ 
stone  and  dry  tuff6  with  similar  density  and  porosity)  and  the  average 
initial  bulk  modulus  (from  Figures  3  and  7)  are  shown  in  Table  3. 


TABLE  3.  INITIAL  ELASTIC  CONSTANTS 


Shear  Velocity 
mV  nsec 

(1  Wz) 

Longitudinal  Velocity 

Calculated  Frost  Ultrasonic  Heasureatents 

Initial 

Bulk  Modulus 
Frost  Hydrostat 
K,  kbars 

MV'wsec 

(1  MHz) 

Constrained 

Modulus 

1,  kbars 

Bulk 
fc/dul us 

K,  ktars 

Young’s 
Modulus 
Y.  kbars 

Poisson 

Ratio 

T 

Shear 
Hodul us 
6,  knars 

Dry  NTS  Tuff 

1.52 

2.63 

117 

64.8 

97 

.25 

39 

36 

Dry  Sandstone 

1.C3 

2.20 

1 

97 

34.6 

97 

.03 

47 

■foot’ 

14 

“linear" 

118 

The  average  slope  of  the  hydrostat  in  the  approximately  elastic 
region  does  not  agree  with  that  obtained  from  ultrasonic  velocity 

g 

measurements.  This  can  be  seen  in  Figure  10,  where  the  slope 
obtained  from  ultrasonic  measurements,  labeled  ultrasonic  K,  has 


Figure  10.  Hydrostatic  response. 


As  the  hydrostatic  pressure  is  increased— beyong  the  approxi¬ 
mately  elastic  region— massive  crushup  of  the  rock  occurs.  The 
rock  has  "yielded"  under  bulk  hydrostatic  pressure.  Crushup 
continues  until  the  porosity  has  been  eliminated.  However,  the 
pressure  volume-change  curve  steepens  considerably  well  before  all 
porosity  is  removed.  It  is  clear  from  Figure  10  that  high  mean  normal 
stresses  are  required  to  eliminate  the  porosity;  i.e..  for  the  hydro- 
stat  to  approach  closely  the  "solid"  curve  (taken  approximately  as 
quartzite)  offset  by  a  volume  strain  equal  to  the  initial  porosity. 

Hydrostatic  Compression  Plus  Shear  Stress 

The  addition  of  shear  stress  to  a  hydrostatic  compression  in 

highly  porous  rocks  has  the  effect  of  Increasing  the  volume  strain  over 
2  8-10 

some  stress  range.  *  That  Is,  stress  can  help  to  collapse  the 
structure,  and  more  porosity  Is  removed  for  a  given  mean  normal  stress. 

This  Is  shown  clearly  in  Figure  11  where  all  load  paths  containing 
shear  stress  lie  on  the  compacted  side  of  the  hydrostat  for  some 
stress  region.  In  this  figure  the  curves  marked  with  o3  refer  to  the 
trlaxial -stress  tests  and  the  curve  marked  to  the  proportional -stress 

tests. 

Dilation  (which  occurs  under  selected  load  paths^)  tends  to  exceed 
compaction  as  the  specimen  nears  the  peak  shear-stress  It  will  support. 

In  the  constant  confining  pressure  tests.  This  dilation  gives  a  volume 
Increase  more  than  offsetting  the  additional  compaction  achieved  under 
compression  at  constant  confining  pressure.  Surprisingly,  dilation  occurs 
when  all  three  stresses  are  compressive,  even  for  materials  that  contain 
large  porosity. 


Shear  Strains  for  Various  Loading  Paths 


The  shear-stress/shear-strain  curves  for  various  loading  paths 
are  shown  in  Figure  12.  An  increase  in  apparent  shear  modulus  and 
shear  strength  with  pressure  is  seen.12*13  Work  hardening  occurs 
over  a  large  strain  range;  however,  it  should  be  noted  that  if  the 
stress-strain  curves.  Figures  4  and  8,  were  replotted  in  terms  of 
true-stress  and  natural -strain,  the  apparent  increase  in  shear 
strength  at  high  shear  strains  would  be  much  less. 


Figure  12.  Shear  strains  for  various  loading  paths. 
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Yield  and  Fracture 


Generally,  rocks  tend  to  behave  in  a  brittle  manner  up  to  some 
confining  pressure  and  then  undergo  a  transition  to  a  ductile  type 
crushup  with  further  increase  in  mean  normal  stress.14  Both  brittle 
and  ductile  need  to  be  qualified,  though.  Brittle  behavior  means  a 
relatively  abrupt  decrease  in  load-carrying  ability,  and  results  in 
a  fractured  sample.  The  stress-strain  curve  generally  does  not  show 
an  instantaneous  drop  of  stress,  but  instead  may  show  a  decrease  in 
load-carrying  ability  with  increasing  strain.  Ductile  behavior 
refers  to  the  ability  to  undergo  appreciable  strain  without  any 
abrupt  drop  in  load,  analogous  to  plastic  flow,  but  generally  not 
"plasticity"  in  the  sense  of  dislocation  motion.  After  ductile  flow 
a  sample  may  be  highly  fractured,  as  for  the  case  of  the  sandstone 
here,  or  it  may  be  intact  with  relatively  little  visual  evidence  of 
crushup,  as  for  the  case  of  the  tuff. 

The  brittle  fracture  and  ductile  crushup  may  be  indicated  in 
stress-space  as  shown  in  Figure  13.  Many  tests  are  required  to  define 
even  approximate  details  of  the  initial  yield  surface  (whether  it  be 
brittle  failure  or  ductile  failure)  ana  its  subsequent  displacement  due 
to  work-hardening.  Some  estimates  can  be  made,  however,  from  Figure  13. 
The  initial  yield  under  hydrostatic  compression  for  this  tuff  has  been 
defined  by  point  A,  point  B  defines  the  yield  point  under  unconfined 
compression,  and  point  C,  the  yield  under  uniaxial -strain  loading. 

These  three  points  approximate  the  initial  yield  surface;  however,  this 
is  at  best  only  as  good  as  the  ability  to  distinguish  the  yield  during 
each  of  these  tests. 


SHEAR  STRESS  (CT,  “  cr3)/2,kbars 


► 


► 


r 


i 


Figure  13.  Brittle  fracture  and  ductile  crushup  for  tuff. 


Less  information  exists  regarding  subsequent  yield  surfaces. 

The  data  show  that  if  hydrostatic  loading  occurs  to  point  D,  unloading 
and  reloading  produce  (at  least  approximately)  elastic  behavior 
until  point  D  is  reached  again  (see  Figures  3  and  7).  Hence  the 
yield  surface  has  been  moved  from  point  A  to  point  D. 

Loading  along  a  path  to  point  E  (Figure  13)  is  well  into  the 
inelastic  region.  Unloading  and  reloading  a  "small"  load  increment 
generally  produces  elastic  behavior  until  point  E  has  been  reached 
again  during  reloading.*®  Unloading  and  reloading  a  "large"  load  in¬ 
crement  (say  all  the  shear  load)  does  not  necessarily  produce  elastic 
behavior.1®  Hence,  the  subsequent  yield  surface  after  reloading  is 
not  well  defined,  and  many  more  tests  than  those  conducted  here  are 
required. 
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Tests  have  shown  a  general  path  independence  of  the  brittle 

c 

failure  shear  strains  of  low-porosity  rocks,  and  some  evidence  exists 
to  suggest  shear  stress  limit  in  thi  ductile  region  (Figures  2  and 
6)  is  also  relatively  path  independent.5’1®  However,  this  path  inde¬ 
pendence  does  not  apply  if  plastic  deformation  occurs  on  one  side  of 
the  brittle-ductile  transition  and  the  maximum  stress  is  reached  on 
the  other  side.1® 

The  path  independence  (or  dependence)  of  strains  has  not  been 
so  well  defined,  however.  Some  data  suggest  only  a  weak  path  de¬ 
pendence  for  loading  just  in  a  range  where  rocks  fail  brittlely.17 
This  is  the  consequence  of  the  small  inelastic  deformation  which 
occi’^s  prior  to  brittle  failure.  The  data  here,  in  the  ductile  ranges, 
suggest  a  relatively  strong  path  dependence  of  strains  due  to  large 
inelastic  deformation.  This  is  shown  when  Figure  14  is  combined  with 
Figure  2.  For  example,  the  4  kbar  constant  confining  pressure  test 
crosses  the  uniaxial-strain  loading  path  at  a  stress  difference  of  4 
kbars  (see  Figure  2).  At  this  point  oj  =  8  kbars,  o2  =  o3  *  4  kbars 
for  both  tests.  From  Figure  14  it  is  apparent  that  the  axial  strains 
are  not  the  same  for  both  tests  at  =8  kbars,  nor  are  the  shear  or 
volume  strains.  The  same  reasoning  can  be  made  by  selecting  similar 
strain  states  from  Figure  14  and  noting  the  different  stress  states  in 
Figures  14  and  2  (using  an  interpolation  of  stresses).  Similar  behavior 
was  observed  for  the  sandstone. 


AXUU.  STRAIN  («, ).  * 


Figure  14.  Path  dependence  for  tuff. 

Dynamic  Effects 

Low-to-intermediate  strain-rate  data  and  plane  shock-wave  data 
exist  for  tuff.18,19  These  data  show  a  7  percent  increase  in 
failure  stress  per  decade  increase  in  strain  rate,  but  are  for 
material  with  an  average  density  of  1.45  gm/cm3:  comparison  with 
the  static  response  of  1.7  gm/cm3  samples  would  therefore  have 
little  meaning.  There  do  exist,  however,  shock-wave  data  on  dry 
Mixed  Company  sandstone  from  the  same  location  as  the  material  used 
in  the  static  tests.18  These  data  are  from  dynamic  plate-impact 
experiments  which  produce  conditions  of  uniaxial  strain.  The  dynamic 
uniaxial- strain  response  from  two  plate-impact  experiments  is  compared 
with  the  static  data  in  Figure  15.  It  can  be  seen  that  there  is 
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considerable  variation  in  the  shock-wave  data:  therefore,  firm  con¬ 
clusions  regarding  the  magnitude  of  dynamic  effects  are  difficult 
to  make. 


VOLUME  STRAIN,  )L,% 


Figure  15.  Comparison  of  SRI  plate-impact  data 

(reference  18  )  and  static  uniaxial  strain 
data  on  dry  Mixed  Company  Sandstone. 
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IV.  A  THEORETICAL  CONSTITUTIVE  MODEL 

Rocks  generally  possess  a  significant  degree  of  porosity  and, 

as  has  been  shown  in  Sections  II  and  III,  exhibit  pore  collapse 

under  hydrostatic  pressure  (a  few  tenths  of  a  kilobar  here).  This 

behavior  is  in  some  ways  similar  to  that  of  porous  ductile  metals 

20-22 

subject  to  hydrostatic  loading.  However,  in  ductile  metals, 

an  increase  in  mean  normal  stress  results  in  pore  compaction  regard¬ 
less  of  the  deviatoric  contribution,  while  rocks  can  undergo 
volne  expansion  in  certain  circumstances  for  which  the  mean  stress 
is  Increasing.  These  effects  of  volume  compaction  and  volume 

dilation  have  been  mathematically  described  with  reasonable  success 

23-26 

for  a  number  of  geological  materials  by  the  so-called  'cap  model." 

These  models  are  all  quite  sophisticated  and  generally  change  in 

form  depending  on  the  material  that  is  being  described. 

The  Important  question  that  needs  to  be  answered  is  what  is 

the  simplest  model  that  can  be  used  to  represent  the  main  features 

of  rock  behavior?  To  attempt  to  answer  this  question  and,  at  the 

same  time,  obtain  constitutive  relations  for  the  two  materials 

Investigated  here,  a  special  case  of  the  cap  model  described  by 

26 

Green  and  Swanson  for  concrete  is  used. 

This  model  can  be  defi’tcd  in  terms  of  a  pair  of  yield  surfaces, 
which  are  most  conveniently  represented  in  terms  of  Jj  and  J2,  the 
first  and  second  Irvarlants  of  the  stress  tensor: 

Jl  ■  On  +  022  +  °33»  (1) 
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J2  =  g-  [(on  -  O22)  +  (o22  “  C  33)  +  (033  “  °ll)  J 

+  a\ 2  +  023  +  o 31 •  (2) 

Figure  16  shows  the  relative  positions  of  the  two  yield  surfaces  in 
the  J i,  /J7 plane  before  (solid  lines)  and  after  (dashed  lines)  a 
tri axial  compression  test  ACD.  These  surfaces  are  defined  in  terms 
of  the  functions  Fj  and  F2  according  to 

F 1  (Jj,  XQ  -  Ki  -  0.  (3) 

F2  (Ji,  -  K2  =  0,  (4) 

where  K,  and  K2  are  coefficients  that  increase  with  the  amount  of 
plastic  flow  and  hardening  that  has  taken  place.  The  initial  values 
of  Kj  and  K'2  are  Ki  and  K2,  respectively.  The  surface  F j  -  K2  =  0 
determines  the  onset  of  plastic  shear  deformation,  and,  for  the  case 
of  positive  slope  as  shown  in  Figure  16,  is  the  source  of  volume 


a  c  a  Ji 


Figure  16.  Initial  yield  surface 

positions  for  cap  model. 
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dilatancy  when  an  associated  flow  rule  is  used  to  define  the  plastic 
strain  components.  Also,  F2  -  K2  =  0  defines  the  cap,  which  is 
necessary  to  account  for  volume  compaction.  For  example,  irreversible 
volume  compaction  occurs  in  hydrostatic  compression  when  point  B  on 

the  axis  in  Figure  16  is  reached. 

Because  the  yield  surface  consists  of  two  independent  segments, 

the  plastic  strain  components,  as  determined  by  an  associated  flow 

26 

rule,  are  given  by  Green  and  Swanson 

,  p  . ,  x  (5) 

dEu  -  x>  x?  • 

where  a  .  is  a  component  of  the  stress  tensor  and  and  \2  are  as 
•  J 

yet  undetermined  multipliers.  The  total  strain  increment  is  given 
by  the  sum  of  dejj  and  an  elastic  strain  obtained  from  Hooke's  law: 

“u  1  <i +  ‘  (1/6|‘  ’ 1/910  6lj  “ Jl’  <6> 

where  K  and  ,  are  the  bulk  and  shear  moduli,  and  sfj  is  the  Kronecker 
del ta. 

Each  of  the  two  yield  surfaces  is  assumed  to  harden  independently 
according  to  the  following  relationships: 

dk,  -  U/hj)  dVif,  (7) 

dk2  =  Cl/ha)  def,c»  ^  ^ 


where 


[(*?,  -  d&>2  -  (*8*  •  d‘P33>2  +  (dEP33  - 


+  (de^)2  +  (d£?3)2  +  ^d£2?)2 
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The  term  de^.  in  Eq.  (8)  is  the  plastic  volume  compaction  obtained 


from  the  surface  F2  -  K2  =  0.  That  is, 

,p  ,  8F2 

deL  *  3X?  tt-  . 


(10) 


VC  3Jj 

The  terms  hj  and  h2  in  Eqs.  (7)  and  (8)  can  be  functions  of  stress, 
plastic  strain,  etc.  In  fact,  it  is  through  the  functions  hj  and  h2 
that  the  de ' ails  of  the  experimental  data  are  fit  by  the  theory. 

It  can  be  shown  with  considerable  algebraic  manipulation  that, 
in  the  region  of  plastic  deformation,  Xj  and  x2  are  given  b^ 


Equation  (11)  holds  for  Fj  -  -  0  and  dFj  >0,  i.e.,  when  the 

stress  state  is  on  the  yield  surface  and  in  a  condition  of  loading. 
Otherwise,  is  equal  to  zero.  Similar  arguments  apply  to  x2  and 
Eq.  (12). 

The  functions  Fj  and  F2  must  be  chosen  in  such  a  way  as  to 
represent  the  essential  features  of  the  experimental  data  without 
unnecessary  complication  in  the  analysis.  This  means  that  we  do  not 
want  to  choose  functions  that  possess  a  great  deal  of  sophistication 
unless  it  is  found  necessary  in  representing  some  important  aspect 
of  the  data  that  cannot  be  described  by  simpler  means.  For  this 
reason  we  choose  the  least  complicated  forms  for  Fj  and  F2  that  will 
allow  both  volume  dilatancy  and  significant  volume  compaction: 


(13) 


Fi  =  "  ®3 1 » 
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h  -  Jl. 


(14) 


Equation  (13)  defines  a  shear  yielding  surface  of  the  Coulomb-Mohr 
type  where  a  is  the  friction  angle  and  Eq.  (14)  defines  a  vertical  cap 
It  will  be  seen  that  these  very  simple  functions  go  a  long  way  in 
representing  the  mechanical  properties  of  the  dry  tuff  (NTS)  and  Mixed 
Company  sandstone  described  in  this  report. 

In  order  to  represent  the  mechanical  parties  of  a  particular 
material,  all  that  remains  to  be  done  is  to  specify  the  elastic 
moduli  and  the  functions  h,  and  h2  in  Eqs.  (7)  and  (8).  This  is  not 
generally  an  easy  task  and  the  functional  fom  of  h,  and  h2  can 
change  from  one  material  to  another.  However,  there  are  a  few 
guidelines  that  can  he  followed  In  choosing  these  functions. 

It  is  advantageous  to  deal  with  the  deviatoric  and  hydrostatic 
effects  separately:  Eg.  (7)  is  related  to  the  deviatoric  response 
end  Eq.  (8)  determines  the  hydrostatic  compaction  properties.  The 
plastic  deviatoric  strain  increment  is  proportional  to  dVIf  in 
Eq-  (7),  which  can  be  written  as 

dV5?=  hj  dFj  =  hj(d^-  a  d  Jj),  (15) 

where  Eq.  (13)  for  the  shear  fail  ire  envelope  has  been  used.  From 
tri axial  tests  on  rocks,  it  is  generally  found  that  the  plastic  strain 
increment  for  a  given  deviatoric  stress  increment  increases  with 
plastic  strain  and  decreases  with  confining  pressure;  c.f..  Figures 

12a  and  12b.  Thus,  these  general  effects  must  be  reflected  in  a  choice 
of  the  function  hj. 
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The  amount  of  irreversible  volume  change  due  to  pore  compaction 
is  controlled  by  Eq.  (8)  which  can  be  written  as 

*;c  =  h2dF2.h2dj1,  ns) 

where  Eq.  (14)  for  the  vertical  cap  has  been  used.  The  function  h2 
is  chosen  to  be  large  where  the  rate  of  plastic  volume  compaction  is 
large,  and  vice  versa.  For  example,  when  the  material  state  approaches 
zero  porosity  in  compaction,  h2  must  vanish.  From  observation  of 
mean  normal  stress  and  volume  strain  in  triaxial  tests,  it  is  observed 
that  the  presence  of  deviatoric  effects  initially  enhances  the  rate 
of  plastic  volume  strain.  All  of  these  phenomena  influence  the  choice 
of  the  function  h2. 

To  obtain  a  fit  to  the  data  on  dry  tuff  (NTS)  and  Mixed  Company 

sandstone  presented  in  this  report,  the  procedure  outlined  in  the 

previous  two  paragraphs  is  followed.  The  actual  choice  of  the  functions 

hj  and  h2,  and  the  numerical  values  which  go  into  the  functions,  is 

made  by  trial  and  error.  This  is  precisely  the  approach  that  was 

26 

used  in  developing  a  constitutive  model  for  concrete.  One  can  start 
with  simplified  functional  forms  for  hj  and  h2  based  on  few  adjustable 
parameters  and  it  very  soon  becomes  apparent  that  these  simple 
expressions  are  not  adequate  for  fitting  all  of  the  data.  Conse¬ 
quently,  the  functional  expressions  become  modified  to  a  considerable 
extent  and  take  on  a  fairly  complex  appearance.  For  the  data  pre¬ 
sented  here,  the  following  functional  forms  for  ht  arid  h2  were  used 
with  independent  variables  Ji,  K},  etc.  in  units  of  kbar: 
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Tuff 


Mkbsr'1) - ^56_exp  [  1#62  K, 

1  +  0.02  J2 

h2(kbar  *)  =  0.0043  J^xpl-  0.075  J, 

4.31  [1  -  exp(-  lOi/JJ)] 

+ - — - - 

1  +  0.038  J ! 

K  =  200  -  164  exp(-  0.071  kbar, 

M  -  52  kbar, 

Kj  =  0,  K2  =  1  kbar, 
a  =  0.064. 

Sandstone 

h ! (kbar- 1 )  =  — — r  exp(2  Ki),  (19) 

1  +  0.1 

h2(kbar-1)  =  0.01  f ( J i ) 

/  £p  \0.3 

+  2 . 75  [l  -  exp (-^37)3  exp [-2.5  f(Ji)]  \1  -  Q^y/  »  (2°) 

f(Ji)  =  0,  Ji  <  10.5  kbar, 

f(Jj)  =  1  -  exp [(10.5  -  Ji)/3],  Ji  -  10.5  kbar, 

K  =  130  -  120  exp(-  Ji/3)  kbar, 

v  *  20  kbar, 

Ki  5^2s  0, 

a  *  0.082. 


/(I  +  0.038  Ji)]  >  (17) 

i*4) 
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Based  on  Eqs.  (17)  and  (18),  and  numerical  integration  of  Eqs. 

(5),  (6),  (11)  and  (12),  calculations  of  the  mechanical  response  of 
tuff  under  hydrostatic  and  triaxial  testing  conditions  are  show,  in 
comparison  with  experimental  data  in  Figures  17  and  18.  A  minor 
shortcoming  to  the  model  for  shear  deformation  of  dr,  tuff  is  the  fact 
that  it  predicts  a  sharp  yield  point  for  the  onset  of  plastic  flow 
and  the  data  (shown  in  Figure  17  )  do  not  indicate  this  behavior.  The 
necessary  codifications  to  the  constitutive  model  in  order  to  include 
the  observed  shear  behavior  are  extensive  and  would  require  completely 
different  hardening  rules.  However,  it  can  be  seen  from  these  figures 
that  the  simple  nodel  based  on  the  linear  Coulomb-Hohr  shear  failure 
envelope  and  the  vertical  cap  does  a  fairly  good  job  in  representing  ell 
the  experimental  data  including  the  phenomenon  of  volume  dilatancy 

as  shown  in  Figure  18. 

Similar  calculations  were  compared  with  experimental  data  on 
Mixed  Company  sandstone  and  the  results  are  shown  in  Figures  19  and 
20  Again,  all  of  the  essential  features  of  the  experimental  data 


have  been  represented. 

In  both  cases  it  is  believed  that  whatever  discrepancies  exist 
between  theory  and  experiment  are  due  to  shortcomings  in  the  functions 
hi  and  h2  rather  than  to  the  choice  of  yield  envelopes.  For  example,  a 
reduction  in  a  at  large  mean  normal  stress  would  result  in  greater 
disagreement  between  theory  and  experiments  in  Figures  18  and  20. 
Because  of  the  rather  cumbersome  functions  necessary  to  fit  the  data, 
it  seems  that  we  have  not  been  judicious  in  the  choice  of  the  parti¬ 
cular  forms.  This  problem  can  be  solved  only  by  additional  investi¬ 
gation  of  other  functions  in  representing  the  mechanical  behavior  of 
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rocks.  Nevertheless,  the  mathematical  analysis  presented  here  seems 
to  be  as  adequate  as  any  work  to  date  In  describing  the  loading  pro¬ 
perties  of  rocks  within  the  limits  of  experimental  observation. 
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Figure  17.  Comparison  of  predicted  and  measured  shear- 
strain  response  for  dry  NTS  ash-fall  tuff. 


Figure  18.  Comparison  of  predicted  and  measured  volume-strain 
response  for  dry  NTS  ash-fall  tuff. 
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Figure  20.  Comparison  of  predicted  and  measured  volume- 

strain  response  for  dry  Mixed  Company  sandstone. 


V.  DISCUSSION  AND  CONCLUSIONS 


The  stress -strain  measurements  obtained  show  that  two  very 
porous  geologic  materials--a  dry  tuff  and  a  dry  sandstone--exhibit 
complex  deformation  and  strength  features.  These  are: 

a.  an  increase  in  shear  strength  with  mean  normal  stress, 

b.  large  permanent  compaction  dependent  upon  both  mean  normal 
stress  and  shear  stress, 

c.  path -dependent  strains, 

d.  strong  strain  strengthening  (work  hardening), 

e.  volume  compaction  is  followed  by  volume  dilatancy  for  loading 
paths  in  which  the  ratio  d(oj  -  a3)/<J(a3)  is  sufficiently 
large  (such  as  the  uniaxial  compression  test  performed  here). 

The  general  features  above  are  exhibited  by  both  rocks.  However, 
observation  of  the  micro-structure  of  each  material  showed  seme  dif¬ 
ferences.  The  tuff  contains  pores  from  Angstrom  size  (10~^  mm)  to 

•  -? 
millimeter  size,  with  much  of  the  porosity  believed  to  be  10  mm  and 

27 

oelow.  Apparently  relatively  little  of  the  total  porosity  exists 
as  open  microcracks.  The  sandstone,  contains  most  of  Its  porosity 
in  the  form  of  holes  about  0.1  mm,  plus  some  amount  of  fracture 
(microcrack)  porosity. 

Elasticity,  even  non-linear  elasticity,  cannot  fit  the  observed 
general  features,  nor  can  simple  plasticity  handle  both  the  compaction 
and  the  dilation  observed.  The  only  theory  that  can  presently  accom¬ 
modate  all  the  above  features  is  the  "cap  model."  This  model  has  been 
fit  to  the  observed  stress-strain  response  using  a  relatively  simple 
hardening  law,  Eqs.  (7)  and  (8).  The  model  handles  the  crushup  under  hy¬ 
drostatic  pressure,  enhanced  compaction  under  various  shear  stress. 
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volume  dilation  under  some  load  paths,  work  hardening  and  the  effect 
of  mean  normal  stress  on  shear  strength.  For  the  hydrostatic  pres¬ 
sure  and  triaxial -stress  load  paths  presented  here,  the  model  fits 
the  stress-strain  response  relatively  well.  To  compare  experiment 
and  theory  for  uniaxial  strain  requires  a  different  calculation  than 
presented  in  this  report.  However,  it  is  expected  that  the  model  pre¬ 
sented  here  should  give  reasonable  agreement  with  the  measured  uniaxial 
strain  response. 

The  most  difficult  features  to  model --regions  that  currently 
need  more  research— are  (i)  the  tensile  and  very  low  mean  normal 
stress  region  and  (ii)  the  hardening  that  occurs  for  various  load 
paths.  A  combination  of  theoretical  modeling  and  experiments  are 
needed  to  extend  the  cap  model  to  adequately  cover  these  two 
regions. 

Additionally,  to  model  in-situ  rock  masses  many  real-world 
in-situ  features  must  be  considered.  Such  real  effects  as  anisotropy 
(tectonic-stress  induced  and  natural)  and  nonhomogeneity,  joints 
and  faults,  in-situ  stress  and  moisture  content  and  pore  pressure 
must  be  considered.  Modeling  these  features  was  not  the  intent  of 
this  work;  some  of  the  effects  can  be  addressed  in  the  laboratory, 
while  others  will  require  field  work  experimentation. 
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DETAILED  EXPERIMENTAL  TECHNIQUES 


Stress  and  Strain  Measuring  Devices 

Figure  A-l  is  a  photograph  of  a  test  assembly  showing  stress  and 
strain  sensors  employed  in  a  typical  test.  An  Important  feature  of 
the  testing  system  described  below  is  that  all  principal  stress  and 
strain  components  of  the  specimen  are  measured.  Measurements  are 
performed  in  a  high  hydrostatic  pressure  environment.  Hydrostatic 
pressure,  a3,  is  monitored  by  a  350  ohm  manganin  wire  coil  connected 
electrically  as  a  single  active  leg  of  a  Wheatstone  bridge.  Stress- 
difference,  (o i  -  a 3) ,  is  monitored  by  a  strain-gaged  load  cell  of 
conventional  design  in  the  form  of  a  high-strength  steel  cylinder. 
Calibration  of  these  transducers  is  performed  using  precision  laboratory 
standards  for  pressure  (a  7  kbar  Heise  gage)  and  force  (a  6400  kilobram 
proving  ring). 

The  strain  transducers  measure  the  axial  strain,  c1(  and  trans¬ 
verse  strain,  e2  *  e3.  The  axial  strain  transducer  which  can  monitor 
a  30%  axial  shortening  strain  (reference  length  is  3.81  cm)  with  an 
accuracy  of  0.16%  strain,  is  shown  in  Figure  A-2.  There  are  four  thin 
cantilevered-beam  arms  which  are  strain-gaged  (350  ohms)  on  both  flat 
sides  and  wired  into  a  four-legged  Wheatstone  bridge  with  two  active 
gages  in  each  leg.  These  four  arms  are  rigidly  attached  to  the  base 
ring  which  is  mounted  to  the  hardened  steel  specimen  lower  endcap. 

When  displacement  occurs  in  the  specimen,  the  displacement  ring 
mounted  to  the  upper  endcap  displaces  the  cantilever  arm  tips  along  a 
conically  tapered  surface.  The  resulting  strain-gage  bridge  unbalance 
is  self-compensating  for  the  small  rotations  of  the  endcaps  sometimes 
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Figure  A-l.  Specimen  transducer  assembly. 
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encountered  while  hydrostatically  loading  a  specimen.  Calibration  over 
the  entire  strain  measurement  range  of  the  transducer  is  performed 
using  a  laboratory  precision  height  gage. 

The  transverse  strain  transducers  (establishing  the  orthogonal 
axes  for  z2  and  £3)  are  cantilevered-beam  systems  also.  Two  pairs 
of  cantilevered  arms  each  carry  two  strain-gages.  Opposite  arms  are 
wired  into  a  single  Wheatstone  bridge  permitting  measurement  of 
diametric  displacement.  Thus  the  transducer  relies  upon  the  spring 
force  of  the  cantilever  arms  to  sense  the  specimen  diameter.  The 
contactor  screws  allow  the  arms  to  be  pre -strained  by  an  amount 
sufficient  to  maintain  contact  during  a  diametric  strain  of  about 
8%  (for  19  mm  reference  diameter).  Calibration  of  these  cantilevers 
is  performed  in  a  manner  similar  to  the  axial  transducer. 

Two  sets  of  strain  transducers  were  developed,  each  for  a 
specific  test  type.  For  hydrostatic  pressure- volume  and  tri axial 
tests  a  set  of  axial  and  transverse  cantilevers  was  used  to  track 
the  large  volume  strain  (up  to  approximately  ±8%  for  each  axis)  with 
an  accuracy  of  .10%  strain. 

A  second  set  of  transverse  strain  transducers  was  developed 
exclusively  for  static  uniaxial -strain  tests.  The  sensor  illustrated 
1  r.  Figure  A-3  was  used  to  monitor  and  control  diametric  strains  at 


Figure  A-3.  Transverse  strain  cantilever  system  used 
for  uni  axial -strain  tests. 
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and  near  zero.  This  transducer  has  a  strain  range  of  ±1%  of  the  19  mm 
reference  diameter  and  can  be  used  to  maintain  zero  strain  with  an 
accuracy  of  .005%  (50  ^strains)  for  homogeneous  materials.  For  large 
porosity,  non -homogeneous  rocks,  this  high  accuracy  was  not  attained. 
Aside  from  the  non-uniformity  of  the  distribution  of  pores  and 
granular  materials  in  the  rock,  the  sample  surface  plays  an  important 
role.  The  urethane  jacketing  material  is  known  to  intrude  into 
openings  with  diameters  of  the  order  of  .1  to  .5  mm  at  high  pressures. 
This  intrusion  necessarily  reduces  the  jacket  thickness  locally, 
thereby  causing  uncertainty  in  the  jacket  effect  on  the  pretest  cali¬ 
bration.  Our  best  estimate  of  the  Measurement  and  control  of  zero 
lateral  strain  is  an  accuracy  of  .05%  strain. 


Small -Strain  Approximation 

Throughout  this  work  we  have  estimated  the  volume  strain  accord¬ 
ing  to  the  small-strain  approximation:  ey  =  z  e^.  The  exact  expres¬ 
sion  for  volume  strain  includes  second-  and  third-order  terms  in  the 
component  strains.  That  is. 


AV/V0 


+  z 
ij 


€i  ej  +  ej  e2  e3- 


(A-i) 


Let  us  define  the  higher  order  terms  in  Eq.  (A-l)  as  E,  the  error 
involved  in  making  the  approximation  aV/Vq  *  Ey.  Values  of  E/ey  are 
then  given  in  Table  A-l  for  the  maximum  pressures  achieved  in  hydro¬ 
static  tests  and  for  maximum  volume  strains  achieved  in  triaxial  tests 
presented  in  this  report.  The  error  (<  7%)  introduced  by  this 
approximation  is  within  the  experimental  scatter  caused  by  sample 
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TABLE  A-l.  ERROR  INTRODUCED  BY  THE  SMALL-STRAIN  APPROXIMATION 


Material 

°i 

(kbar) 

03 

(kbar) 

E/ev 

(S) 

Tuff 

8.0* 

8.0* 

7.0 

II 

5.0 

2.0 

2.5 

li 

7.8 

4.0 

4.2 

II 

16.3 

8.0 

7.0 

Sandstone 

8.0* 

8.0* 

5.0 

II 

3.7 

1.9 

1.0 

II 

8.5 

4.0 

5.0 

*  hydros tat 

variability.  However,  larger  errors  result  when  the  small  strain 
approximation  is  used  in  regions  of  high  volume  dilatancy.  This  has 
negligible  influence  on  the  experimental  and  theoretical  results 
presented  here. 

Transducer  Pressure  Effects 

When  strain-gages  are  exposed  to  high  pressures,  a  component 

of  the  strain  measurement  is  a  pressure-induced  apparent  strain.  This 

28  29 

pressure  effect  on  single  gages  has  been  described  in  the  literature. 

Tests  have  shown  that  the  pressure-induced  effects  are  independent  of 

high-pressure  electrical  feed-throughs,  in-vessel  lead  connectors, 

28 

solder  connections  and  external  circuitry.  Also,  apparent  strain 
due  to  pressure  was  found  to  be  independent  of  pre-set  strain  condi¬ 
tion.  Tests  performed  in  persuance  of  a  solution  to  this  problem 
have  shown  that  the  repeatability  of  the  apparent  strain  output  for 
a  single  gage  in  the  high  pressure  environment  depends  on  two  para¬ 
meters:  (i)  the  thickness  of  the  epoxy  bonding  layer  and  (ii)  the  rate 
of  pressure  application.  A  method  has  been  developed  to  reproducibly 
attach  strain-gages  with  a  bonding  thickness  of  '3  y.  Tests  are  now 
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performed  at  the  pressurization  rate  of  about  1/4  kbar/minute  for 
hydrostatic  tests.  Particular  care  is  paid  in  placement  of  the 
strain  gages. 

Ten  Ki lobar  Testing  System 

A  block  diagram  of  the  10  kbar  system  used  for  the  tests  des¬ 
cribed  here  is  shown  in  Figure  A-4.  The  two  main  press  components 
are  the  loading  system  and  the  pressure  intensifier  system.  Figure 
A-5  is  a  photograph  of  the  loading  frame  showing  the  reaction  platens, 
support  columns,  tie  rods  and  hydraulic  loading  actuator.  The  test 
vessel  is  suspended  from  the  upper  platen  by  a  split  support-ring. 

A  one-inch  diameter  loading-piston  is  attached  to  a  104,000  kg  actuator 
and  is  guided  by  the  loading-piston  seal  surface  of  the  test  vessel. 


Figure  A-4.  Block  diagram  of  10  kbar 
testing  machine. 


30- 
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A  5-cm  diameter  base  plug  supports  the  specimen,  strain  sensors, 
pressure  coil,  load  cell  and  contains  a  Bridgman-type  seal  and  24  high- 
pressure  conical  electrical  feedthroughs.  Figure  A-l  shows  the  base 
plug  with  specimen  sensors  mounted.  The  base  plug  also  transmits 
the  specimen  axial  loading  thrust  through  its  support  spacers  to  the 
platens.  A  small  hydraulic  handjack  is  used  to  overcome  the  base  plug 
seal  friction  and  insert  the  base  plug  test  assembly. 

The  pressure  intensifier  ram  is  attached  to  a  107,000  kg  hydraulic 
actuator  and  with  a  single  full  stroke  produces  a  pressure  In  excess 
of  11  kbars  in  the  test  system.  The  nower  plug  contains  a  moving 
high-pressure  seal  and  rides  freely  on  the  intensifier  ram.  For  ease 


Figure  A-5.  10  kbar  system  loading  frame. 


of  operation  a  low  pressure  pump  may  also  be  used  to  pre-charge  the 
system  (up  to  about  0.7  kbars)  as  well  as  fill  the  test  vessel.  Both 
the  test  vessel  and  the  pressure  intensifier  vessels  are  4340  steel 
that  have  been  autofretteged  (strain  strengthened)  to  a  condition  that 
prevents  further  yielding  at  10  kbars  internal  pressure. 


Servocontrolled  System 

Both  press  components  are  servocontrolled  as  illustrated  in 
Figure  A-6.  The  details  of  selected  type  of  sensor  feedback,  command 
function,  loading-stress -rate  or  strain-rate,  etc.,  are  determined  by 
the  type  of  test  performed.  Categories  of  tests  include  (a)  hydro¬ 
static,  (b)  uni  axial -strain,  (c)  confined  triaxial -compression  or 
other  stress/strair.  histories. 
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a.  Hydrostatic  tests  are  performed  in  which  the  three  principal 
strains  are  recorded  plus  the  fluid  pressure  (symnetrical  stress 
loading).  For  these  tests  the  command  function  generator  produces  a 
symnetrical  sawtoothed  voltage  ramp  (frequency  less  or  equal  to  0.001 
Hz)  for  the  intensifier  servocontroller.  The  pressure  feedback  monitor 
conditions  the  servo  output  signal  for  a  null  offset.  Strains  and 
pressure  data  are  recorded  during  the  ramp  cycle. 

b.  Uni  axial -strain  tests  are  performed  in  the  following  manner: 
a  pretest  run  (no  test  specimen  within  the  vessel)  to  maximum  pressure 
is  performed  while  recording  the  pressure-effects  of  the  transverse 
strain  sensors  on  an  analog  chart  recorder.  The  command  function  gener¬ 
ator  is  set  to  provide  the  voltage  ramp  (frequency  less  or  equal  to 
0.001  Hz)  for  either  the  intensifier  or  load  controller.  If  the 
pressure  servocontroller  is  selected,  the  load  actuator  is  manually 
controlled  and  vice  versa.  As  the  pressure  increases,  load  is  applied 
to  the  specimen  to  just  compensate  for  the  pressure-induced  compres¬ 
sive  diametric  strains.  When  zero  transverse  strain  is  achieved, 

the  transverse  strains  sensors  output  follows  exactly  along  the  pretest 
tracing. 

c.  In  triaxial -compression  tests,  a  hydrostatic  loading  test 
is  initially  performed  to  the  desired  confining  pressure.  The  command 
function  generator  is  then  set  to  provide  voltage  ramp  (frequency  less 
or  equal  to  0.001  Hz)  for  load.  As  load  increases,  the  three  principal 
strains,  fluid  pressure  and  the  load  (stress  difference)  are  recorded. 
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Data  Acquisition  System 

figure  A-7  shows  the  data  acqoisitioh  system  used  for  the  tests 
described  below.  A  Digital  Equipment  Corporation  PDP  Lab  11/20 
computer  is  the  primary  recording  device  and  provides  real-time  test 
data  monitoring  and  analysis.  Each  test  sensor  provides  input  to  one 
channel  of  an  8-channel,  12-bit,  analog-to-digital  converter  which 
operates  over  a  10  volt  range  with  a  maximum  resolution  of  2.4  milli¬ 
volts  (equivalent  to  5  microvolts  out-of -bal ahce  voltage  at  the  bridge). 
The  computer  console  contains  the  central  processor,  core  memory,  a 
two-color  CRT  plotting  screen,  a  high-speed  paper  tape  punch  and  read 
aid  the  A/D  converter.  Peripherals  are  addressable  in  real-time  which 
allows  instantaneous  review  of  acquired  data  and  test  progress.  Both 
the  CRT  plotter  and  the  Valtech  digital  plotter  are  used  to  display 


Figure  A-7.  Data  acquisition  system. 
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The  instrumentation  console  shown  in  Figure  A-8  contains 
Wheatstone  bridge  current  supplies,  hookup  patchboards,  balancing 
networks,  calibration  shunt  networks  and  output  amplifiers  for 
eight  channels  of  data.  Normally  only  five  channels  are  used; 
pressure,  load,  axial  strain  and  two  transverse  strains.  Data 
monitoring  using  X-Y  recorders  as  output  terminals  is  also  available 
at  the  testing  machine  servo-control  console. 


Figure  A-8.  Instrumentation  console. 
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DETAILED  MATERIAL  DESCRIPTION  AND  SPECIMEN  PREPARATION 


The  Nevada  Test  Site  tuff  samples  were  obtained  from  the  left 
rib  of  tunnel  U12e.06  between  stations  #9  +95  and  #10  +50.  This 
tunnel  is  located  in  the  Rainier  Mesa,  Nevada  Test  Site,  Mercury, 
Nevada,  within  ash-fall  tuff,  containing  both  reworked  ash-fall  and 
tuffaceous  sandstone.  The  samples  are  mottled  pink  and  are  mainly 
fine  grained  with  grains  of  less  than  1  mm,  but  contain  occasional 
coarser  fragments  of  pumice  granules,  lithic  grains  and  phenocrysts 
greater  than  2  to  5  mm  in  size.  Figure  B-l  shows  a  photomicrograph 
at  50  times  magnification.  The  glass  shards  that  originally  formed 
most  of  the  matrix  were  crystallized  as  zeolites,  cristobalite  and 
possibly  clay.  Under  standard  room  conditions,  the  density  for  the 


P- Pores 

PF- Pumice  Fragments 
Q-Quartz 


Figure  B-l.  Polarized  reflected-light  photomicrograph 
of  ash-fall  tuff  showing  intragranular 
porosity  and  microporosity  inside  pumice 
fragments . 
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samples  is  approximately  1.7  gm/cm3  with  grain  density  of  approximately 
2.46  gm/cm3.  Porosity,  as  derived  from  the  relation 


n  =  ( 1  -  p  d/p  g)  x  100 


where  pd,  Pg  are  the  dry  and  grain  density  respectively,  is  31  percent. 

The  other  porous  rock  tested  in  this  program  is  from  the  Mixed 
Company  test  site  near  Grand  Junction,  Colorado.  It  is  a  massive, 
medium-grained,  buff-colored  sandstone.  Microscopic  examination 
indicates  it  is  composed  of  almost  100%  rounded  quartz  grains  bonded 
with  silica.  Figure  B-2  shows  a  photomicrograph  at  50  times  magnifi- 
cation.  Grain  size  is  0.5  to  1.0  mm.  Variations  in  density  (1.94  to 
2.03  gm/cc)  can  be  attriboted  to  variations  in  porosity  rather  than 
any  difference  in  mineralogy,  the  rock  is  highly  permeable  and  is 

f 


Figure  B-2. 


Thin  section  photomicrograph  of  Mixed  Company 
sandstone  showing  intragranular  porosity  and 
quartz  grains.  9  stands  for  quartz,  R  for 
rock  fragment  (microcrystalline),  h  for 
pores,  F  for  feldspar. 
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friable  when  wet.  Intact  cores  a  few  feet  in  length  can  be  extracted 
Indicating  the  rock's  massive  nature.  However,  faint  relict  bedding 
planes  are  observed  on  a  macro-scale  within  the  unit. 

Test  specimens  were  obtained  from  a  36-cm  diameter  drill  hole 
in  Area  B  (the  Physics  International  Site)  some  122  meters  south, 
southeast  of  the  working  point  of  a  453. 5  metric  ton  surface  high 
explosive  test.  Typical  dry  laboratory  density  was  2.0  gm/cm3;  the 
corresponding  grain  density  was  2.66  gm/cm3  giving  an  average  porosity 
of  approximately  25%. 

Specimens  used  for  the  tests  were  dry,  solid,  right-circular 
cylinders  1.9-cm  diameter  by  3.8-cm  in  length.  A  1.9-cm  core 
drill  was  used  to  core  the  specimens  from  larger  blocks.  The  ends 
were  ground  with  a  precision  surface  grinder  to  be  as  smooth  as  the 
rock  granularity  allowed  ('.003  cm  rms)  and  parallel  within  +.002  cm. 
The  density  of  each  specimen  was  determined  by  weighing  the  specimen 
(±.01  gm)  and  dividing  by  the  calculated  volume  found  by  averaging 
measurements  of  the  sample  dimensions.  Lengths  and  diameters  were 
measured  accurately  to  ±.003  cm.  By  this  procedure  the  estimated 
density  accuracy  was  ±.01  gm/cm3. 


